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Abstract 



The new particle Z(3930) found by the Belle and BaBar Collaborations through 
\Q ■ the 77 — > DD process is identified to be the Xc2(2P) state. Since the mass of this 

particle is above the DD^ threshold, the OZI-allowed two-body strong decays are 
r*j ' the main decay modes. In this paper, these strong decay modes are studied with 

Q-c two methods. One is the instantaneous Bethe-Salpeter method within Mandelstam 

Q_|, formalism. The other is the combination of the 3 Po model and the former formalism. 

The total decay widths are 26.3 and 27.3 MeV for the methods with or without the 
3 -Po vertex, respectively. The ratio of over r^^, which changes along with the 
mass of the initial meson is also presented. 

■ 1 Introduction 

It is well known that the quark potential models predict abundant charmonium spectra. 
^ | In these spectra, only the radial excitations of the S'-wave spin-singlet and triplet states 
have been observed experimentally. For others such as the four P-wave states 3 Pj and 
1 Pi, only the ground states have been observed and studied carefully. The identification 
of these P-wave radial excitations will greatly support the correctness of potential models 
and supply more playgrounds for the study of QCD non-perturbative properties. 

Recently, the Belle Collaboration [1] found a new particle in the process 77 — > DD with 
statistical significance of 5.3a. The observed mass and decay width are M = 3929 ± 5 ± 2 
MeV and T tota i = 29 ± 10 ± 2 MeV, respectively. Later through the same process while 
with a little higher statistical significance 5.8<r, the Babar Collaboration [2] also observed 
this particle with compatible mass and decay width: M = 3926.7 ± 2.7 ± 1.1 MeV and 
Ftotai = 21.3 ± 6.8 ± 3.6 MeV. The experimental data favor the 2 ++ assignment, of which 
the (helicity) angular distribution of the decay products has a form sin 4 9. However, as 
pointed by Chen et al. [3], the assumption that two P-wave higher charmonia participate 
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as the intermediate state in the process is also compatible with the experimental data of 
two groups. So a careful study of the decay properties of the x c2 (2P) state is still needed, 
which is helpful for the final confirmation of this state. 

Since the mass of this particle is above the DD^ threshold, but under the D*D* 
threshold, the former channels are allowed by the Okubo-Zweig-Iizuka (OZI) rule and 
can be used to estimate the total decay width of this state. As is well known, these 
open-flavor strong decays closely relate to the non-perturbative properties, of which our 
knowledge is rather poor. A complete solution to this problem needs deep understanding 
of the QCD vacuum. Although we expect the lattice QCD calculations will provide firm 
theoretical predictions in the future, now we are forced to construct phenomenological 
models, e.g. the 3 Po model [4, 5, 6]. Possible microscopic models include the flux-tube 
mode [7], Cornell model [8, 9] with a vector confinement interaction, the model in Ref. [10] 
with a scalar confinement interaction, Field Correlator Method [11] et al. 

Here we try to apply the Bethe-Salpeter (BS) method [12, 13] to the two-body strong 
decay processes within Mandelstam formalism [14]. It is well known that this method 
provides a relativistic description for the two-body bound state. Its instantaneous version 
is extensively used to study the properties of heavy mesons [15, 16, 17]. One notices that 
this method has been mentioned in Ref. [18] and used to study light meson decays in 
Ref. [19], where in addition to the pure quark loop diagram, an instanton-induced vertex 
is considered for the processes with vanishing total angular momentum. Here for the 
charmonium case, we will assume the pure quark loop diagram gives the main contri- 
bution. Non-perturbative effects are included in the wave functions of initial and final 
mesons. These wave functions are obtained by solving corresponding coupled instanta- 
neous BS equations. As expected, the relativistic correction is covered, which is especially 
important for the orbital-excited particles. 

In this method we make two approximations. First, we will not consider the inter- 
actions of final mesons, which can be realized by interchanging pions [20]. This maybe 
important for the threshold-nearby decay processes. It is not clear now how to incorpo- 
rate these effects into the BS method. Second, we will not consider the coupled-channel 
effects [21, 22, 23], which may move the pole mass to the physical one. Here we adjust 
the parameter in the potential to get the correct mass value and obtain the wave function 
for the initial meson. 

For comparison, we will also consider the 3 Po model, which assumes that a qq pair 
is created from the vacuum with a quantum number ++ . This model is extensively 
applied to study the strong decays of light mesons [10], heavy-light mesons [24] and heavy 
quarkonia [25]. It is a non-relativistic model, where the simple harmonic oscillator (SHO) 
wave functions are used both for the initial meson in its cm. system and the final particles 
in their moving systems. As mentioned in Ref. [11], this assumption is valid only for small 
relative velocity near the threshold. For P-wave charmonia, more covariant formalisms 
and realistic wave functions are needed if one wants to get a more reasonable result. In 
Ref. [26], we have combined the BS method with the 3 Pq model to make it applicable 
for the OZI-allowed two-body decay processes. There we also proved that this amplitude 
will reduce to the usual form when the non-relativistic approximation is made. The same 
method will be considered in this paper. 

The paper is organized as follows. In Section 2, we present the calculations by using the 
instantaneous BS method within Mandelstam formalism. The pure quark loop amplitude 
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is considered. In Section 3, within the same formalism, we add a phenomenological 3 Po 
vertex to describe the decay mechanism. The results and discussions are given in Section 
4. 



2 BS Method 

The BS equation is written as 

ST 1 (Pi)x P (q)S2 1 (-P2) = i J ^0- 4 v(p ]q ,k) Xp (k), (i) 

where X P (o) is the BS wave function; Si(pi) and S 2 (—p 2 ) are the quark and anti-quark 
propagators; V(P; q, k) is the interaction kernel; p\ and p 2 are the momenta of quark 
and anti-quark, respectively, which are related to the meson momentum P and relative 
momentum q by 

777 ■ 

P^ = —P + (-lY +1 q, for i = 1, 2. (2) 

mi + m 2 

With the definition p ip = and pf ± = pf — ; ^P M , we can write Si(Jpi) as 

A + A~ 
- USiiJpi) = + — s (3) 

p ip - a;, + ze p ip + cjj - te 



where we have used 



(4) 



^ = \im\-p\^ 



and J = 1 for the quark (i — 1) and J = —1 for the anti-quark (z = 2). 

The 3-dimensional form of the BS wave function is defined as <f(q±) = i J -^Xpio}- 
With instantaneous approximation, V(P; q, k) ~ V(P; q±, k±), Eq. (1) can be written as 

X P (q) = Si(pi)ri P (q ± )S 2 (-p2), (5) 

where 

/d 3 k 
-j—^V(q ±J k ± )<p p (k±). (6) 

Within Mandelstam formalism, the leading order transition amplitude for the strong 
decay (see Fig. 1) can be written as 

(P 1 P 2 \S\P)bs = (2tt) V(P -P 1 - P 2 )M bs 

(2tt) 4 (2-kY ( 27 r) 2 4 Tr[Xf (g)lS2 " 1( ~ P2)(27r)4<54(p2 ~P 22 )^fe) 
x 5r 1 (j9 21 )(2vr) 4 5 4 (p 21 +Pi 2 )Xp 1 (gi)5r 1 (pi)(2 7 r) 4 ( 5 4 (p 1 - Pll )] (7) 

= (2 7 r) 4 5 4 (P-P 1 -P 2 )C / | ^jltlXpCgJ^H-ft^Cfe)^ 1 ^ 

x XPi^O-ST^pi)], 
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Figure 1: Feynman diagram of OZI-allowed two-body decay process 
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Figure 2: Feynman diagram of OZI-allowed two-body decay process with a 3 P vertex. 

where Cf = is the color factor and Xp(<z) = 7°xt(<?)7 - The meaning of is explicit 
from the diagram and they are related to the meson momentum P, and relative momentum 
qi by the same relation as in Eq. (2). The relative momenta of final mesons are related 
to that of the initial meson by = q + (— l) i+1 (ctiP — ct^Pj), where 

rrti ma 

on = ; , au = ■ , for z = l,2. (8) 

mi + m 2 ma + m i2 

From the Eq. (7), we get the Feynman amplitude, 

M BS = C f J ^Trtx^g)^-^-^)^^)^- 1 ^)^^?!)^ 1 ^)] 
f d 4 q 

= C f J j2^ TT l s i(Pi)VpS2(-P22)v P2 S 2 ( y -pi2)fi Pi } 
„ - r [ d 4 q At(q±) At(q±) - , Ajifex) 



(27r) 4 pi p -u;i + ie p p 2 p-^2 + «e P2 £>2i P - ^21 + «e 

Pl2 P - W12 4 

d 3 <f 



(9) 



In the second equation, we have used Eq. (5), and in the third equation, we have inserted 
Eq. (3). For simplicity, we only keep the positive energy projectors which give the main 
contributions. In the parentheses we have used A± 2 (qi±.) = Aj 1 (g 2 ±)- Considering the BS 
equation [15] 

(M — U\ — cu 2 )v + p + (q ± ) = Af Vp (q ± )At, (10) 

(y? ++ is defined as A^-^^-pA^ ) and by doing the residual integral, we get the fourth 
equation. In the fifth equation we have used Eq. (6). 

The wave functions of 2 + , 1~ and 0~ states have the following forms [17, 15, 16], 



^= M ^ + 9^9 3 % + 9 S^h\ (11) 



q 2± ■ e)h 7 + -^{t-ziii i _ ^2?2± • e)h 



where fi, gi and ^ are functions of q, q\ and q 2 , respectively. 

To perform the integral in Eq. (9), we have chosen P 1 as the reference direction, 

Pi -<f =1^11^ cos l5 P 1 -k= | cos 6" 2 , 
q-k=\q\\k\ cos 9, (12) 
cos 9 = sin 8i sin 9 2 cos(0i — <p 2 ) + cos 9i cos 9 2 . 

By defining P^ ± = — -^gjr-P^ and = g^ — , we can express the integrals 
which have free Lorentz indexes as follows, 

J d 3 qd 3 kq» ± F(q ± , k ± ) = f u P^, J d 3 qd 3 kklF(q^ k ± ) = g n P? ± , 

J d 3 qd 3 kq^ ± F(q ± , k ± ) = f 21 P^ ± + fag?, 

J d 3 qd 3 kkW ± F(q ± , k ± ) = g 2l P^P^ + g 22 gT, (13) 
J d 3 qd 3 kq^ ± q a ± F(q ± , k ± ) = /ai^xlAx + f^P^T + P?±9? + A±<7f), 
J d 3 qd 3 ke ± q» ± q a ± F(q ±J k ± ) = g^P^P^P^ + g^P^T + <?3 3 (A±<?f + A±<?f ), 



where and gfy are integrals with no free Lorentz indexes, which can be done numerically. 

Inserting Eq. (11) into Eq. (9) and finishing the trace, we can get the transition 
amplitude with the following forms, 

M( X c2(2P) -+ DD) = e^PfPfc, 
M( X c2(2P)^DD*)^e^ 6 e flv e a P^PpP ls t 2 , 
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where t\ and t 2 are form factors, which are the functions of and gij in Eq. (13). These 
amplitudes actually can be constructed by using the momenta and polarization tensor 
(vector) of initial and final mesons. With consideration of parity properties, we can see 
there is a totally antisymmetric tensor in the second amplitude while the first one does 
not include it. 

The decay width is 

r =j&El*l'. (i5) 

A 

where |A| = y/[M - (Mi - M 2 ) 2 ][M - (Mi + M 2 ) 2 ]/2M is the final meson momentum. 



3 Extended 3 P Model 

The usual 3 Po model is a non-relativistic model with a transition operator g J dxipip\ noQl . c i [27]. 
To combine it with BS wave functions, we extend it to the relativistic form —ig J d A xipip 
(a similar form of interaction is also used in Refs. [11, 22]), where g is parameterized as 
2m q r ). m q is the constitute quark (u, d, s) mass. The interaction strength 7 is dimension- 
less, which is roughly flavor independent in the light meson decay processes [25]. Just as 
Ref. [25] did, we will assume this also applies in the heavy meson case. 

Here for simplicity we will not get the value of 7 by fitting decay widths of other 
channels, instead, we take 7 = 0.35, which is the best-fit value for the usual 3 Pq model [27]. 
In this case, we can only get a roughly estimation for our extended 3 Po model. But 
considering that this model will reduce to the usual one if non-relativistic approximation 
is made, we can argue that by using the same interaction strength we get qualitatively 
how large change can be brought by using a covariant formalism and BS wave functions. 
A more reliable quantity is the ratio of partial decay widths of different channels. 

Within Mandelstam formalism we can write the transition amplitude of the OZI- 
allowed open flavor decay process (see Fig. 2) as 

(PxP 2 \S\P) 3Po = (2tt)V(P -P x - P 2 )Ms Po 

= ~ W j (2tt^ (2t0 4 (2j) 2 4 Tr ^^^ 1 ( p2 ^ 27r ) 4(54 ^ 2 -P2 2 )Xp> 2 )x gi 

xSi\p 1 )(2nY5 4 (p 1 -p 11 )] (16) 

= -ig{2*Y5\P - P, - P 2 ) J -^Trtx^g)^^-^)^^)^^)^ 1 ^)] 

From the above equation, we get the Feynman amplitude [26] 
f d 4 q 

Ms Po = -ig J -—^Tr\x P (<l) S 2 (-P2)x P2 (Q2)x Pl {l^ S i (p01 



(17) 
M — uji — UJ2, 



2uj 12 

where g« = <?+(— l)* +1 (ajP— otuPi). To get the second equation, we have made the residual 
integral and used Eq. (10). The second term in the last parentheses can be neglected for 
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Table 1: The value of — Vo (GeV) for different states. 



States Xc2(2P) 


D+ 


D° 


Df 


JJ*+ JJ*0 


-V 0.177 - 0.098 


0.375 


0.375 


0.432 


0.200 0.200 


Table 2: OZI-allowed two-body strong decay widths (MeV) of x c2 (2P) 
where we have used M = 3930 MeV. 


with two methods, 


Mode D + D- 


D°D° 


DD 


D+D*- 


D°D*° DD* 


Xc2 (2P) (BS) 8.57 


9.47 


18.0 


1.84 


2.33 8.34 


Xc2 (2P) (BS- 3 P ) 10.1 


10.3 


20.4 


1.48 


1.97 6.90 



M Rj uj\ + a; 2 (for large this approximation is not valid, while at the same time the 
wave function is strongly suppressed). Inserting Eq. (11) into above equation we get the 
same forms of the decay amplitudes as Eq. (14). 



4 Results and Discussions 

The potential in Eq. (9) is the Cornell potential. It is the same one as we used to solve 
the instantaneous BS equations. The parameters in the potential are fixed by fitting the 
mass spectra. So we do not introduce new parameters when calculate the decay width. 
In the momentum space, it has the form 

V(q) = (2*) 3 V s (q) + 7 o ® 1° (2ir) 3 V v (q) , 
V s (q) = -(- + Vo)6 3 (q) 



a v " 7r 2 (q 2 + a 2 ) 2 ' 

v[q> ~ 3n 2 q 2 + a 2 -' 

127T 1 



27 ln(a + 



where the following parameter values are used: a = e = 2.7183, a = 0.06 GeV, A = 0.21 
GeV 2 , m c = 1.62 GeV, m s = 0.5 GeV, m u = 0.305 GeV, m d = 0.311 GeV, A QCD = 
0.27 GeV. The value of — Vo is listed in Table 1. For x c2 (2P), we vary it in the range: 
0.177~0.098 GeV. If we let the mass of Xc 2 (lP) coincident with the experimental result, 
we get -V = 0.11 GeV and M(x c2 (2P)) = 3972 MeV, which is about 40 MeV larger than 
the experimental data. To get M(x c2 (2P)) = 3930 MeV, we scale -V to 0.165 GeV. 

The patial decay widths for x c2 (2P) (M = 3930 MeV) with two methods are given in 
Table 2, where DD represents D + D-+D°D° and DD* represents D + D*~ + D°D*° + c.c. 
One notices that the DD* channel contributes about half (a third) of that of the DD chan- 
nel for the BS (BS- 3 Po) model. A precise measurement by future experiments for the ratio 
of two partial widths is expected. The ratio Pr(x c2 (2P) -> P+P")/Pr(x c2 (2P) -> D°D°) 
is 0.90 (for BS), 0.98 (for BS- 3 Po), which is consistent with the experimental value [1]: 0.74 
± 0.43 (stat) ± 0.16 (syst). For the DD* case, Pr(x c2 (2P) -» P + P*")/Pr(x c2 (2P) ->■ 
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Figure 3: (a) IVd, (b) T d q, and (c) T D ^/Y D ^, change along with the mass of the initial 
meson x c2 (2P) 



Table 3: Open-flavor strong decay widths (MeV) of Xc2(2P) state. The second and the 
third columns are our results with BS method and BS- 3 Po model, respectively, where 
the values outside the parentheses are gotten with M = 3930 MeV, while the results in 
the parentheses are calculated with M = 3972 MeV. The fourth column is the usual 3 Pq 
model and in the fifth column, C 3 represents the Cornell coupled- channel model. 



Mode BS BS- 3 P 3 P [27]* C 3 [29] Ref. [3] Exp [28 



DD 18.0 (24.8) 20.4 (24.6) 26 (32) 21.5 

D S D S - (0.347) - (0.620) - (0.5) 

DD* 8.34 (34.7) 6.90 (21.1) 9 (28) 7.1 

total 26.3 (59.8) 27.3 (46.3) 35 (60.5) 28.6 11.98 24 ± 6 

* In Ref. [25], the authors used 7 = 0.4 for the 3 P model with M = 3972 MeV. There they got r tota; 
= 80 MeV. 
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D°D*°) = 0.79 is smaller than that of DD. This is because the former channel has smaller 
phase space, a small mass difference M(D + D*~) — M(D°D*°) will cause a large difference 
of the partial decay widths. 

In Table 3, we present results with different models and experimental data. Our two 
methods get quite close total decay widths, which are also consistent with the experimental 
value. For partial decay widths, BS- 3 Po model gives a larger T D jj, while for T D ^ t , BS 
model gets the larger results. By using the usual 3 P model, Ref. [27] get a large result 
compared with our BS- 3 Po model. As mentioned in Section 3, for we have adopted the 
same interaction strength, the discrepancy may come from using different formalisms 
and wave functions. In Ref. [29], the Cornell coupled- channel (C 3 ) model is used, which 
assumes a current-current interaction with a vector confinement. By using gaussian type 
wave functions and solving the coupled-channel equation, the authors get the complex 
eigenvalue whose real part is the physical mass of the meson and the imaginary part 
related to the decay width. One notices this model gives results very close to ours. 

In Fig. 3, we give the partial decay widths for DD and DD* channels and the ratio 
of the two widths which changes along with the initial meson mass. One can see for 
r^^i, the BS- 3 Po model gives larger results and the difference of two models becomes 
less if the initial meson mass increases. For r^,, the BS model gives larger values and 
the difference of two models becomes larger along with the increasing of x C 2(2P) mass. 
The BS- 3 Po model gives a larger ratio of two partial widths, which decreases about 4 
times when the mass of the Xc2(2P) state increases from 3920 MeV to 3980 MeV. One 
notices that the partial widths with the 3 Po model are 7-dependent, while the ratio is 
7-independent and more reliable. 

By using the widths in Table 3, we can calculate the branching ratio T D jj/T tota i = 
0.684 (BS) and 0.747 (BS- 3 P ). In Ref. [17], the same method gives = 0.534 keV 
for Xc2(2P). Combining the two results, we get r 77 Pr[x c2 (3930) — > DD] = 0.365 keV 
(BS) and 0.399 keV (BS- 3 P ). The experimental data are [1, 2]: 0.18 ± 0.05 ± 0.03 keV 
(Belle) and 0.24 ± 0.05 ± 0.04 keV (BaBar). The theoretical results are consistent with 
the experiment values considering the large uncertainty of the QCD corrections to the 
diphoton decay processes [30]. 

In conclusion, we have used two methods to calculate the OZI-allowed two-body strong 
decay processes of the Xc2(2P) state: the BS method and the extended 3 Po model. In the 
former we did not introduce any new parameters, while in the later a flavor-independent 
interaction strength is used. The total decay width estimated is 26.3 (27.3) MeV for the 
former (later) model with M = 3930 MeV, which is consistent with the experimental 
data. However, more experiments, especially the precise measurements of the partial 
decay widths of DD and DD* channels are still needed for the final confirmation of this 
paritcle. 
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